Short-Time Loschmidt Gap in Dynamical Systems with Critical Chaos 
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We study the Loschmidt echo F{t) for a class of dynamical systems showing critical chaos. Using 
a kicked rotor with singular potential as a prototype model, we found that the classical echo shows 
a gap (initial drop) 1 — Fg where Fg scales as Fg{a,e,i]) — /ci(Xci = ''f~°' oi is the order of 
singularity of the potential, the spread of the initial phase space density and e is the perturbation 
strength. Instead, the quantum echo gap is insensitive to a, described by a scaling law Fg — fq(xq = 
rf /e) which can be captured by a Random Matrix Theory modeling of critical systems. We trace 
this quantum-classical discrepancy to strong diffraction effects that dominate the dynamics. 
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The study of systems with a phase transition was al- 
ways a fruitful subject of study for many areas of theo- 
retical and experimental physics. Specifically in the field 
of disordered metals, the celebrated Anderson Metal- 
Insulator Transition (MIT) [J has been an exciting sub- 
ject of research for more than fifty years. On the other 
hand, the field of quantum chaos brought up a connection 
between quantized chaotic systems and localization ideas 
emerging from solid-state physics It has been shown 
that quantum suppression of classical diffusion is a result 
of wave interference phenomena of similar nature as the 
ones responsible for Anderson localization in disordered 
metals. Quite recently the connection between the two 
fields was further strengthen with the observation that 
certain non-KAM dynamical systems exhibiting classi- 
cally anomalous diffusion, can have statistical properties 
resembling the ones of disordered metals at MIT [^, 0] . 
This phenomenon is referred to as critical chaos. Some 
of these properties include, quantum anomalous diffusion 
[^, multifractal wavefunctions Q, and critical spectral 
statistics Q. Many of these intriguing statistical prop- 
erties can be exactly derived using non-conventional en- 
sembles of random matrices with variance decaying from 
the diagonal in a power-law fashion which in turn 
model a large variety of physical systems [§]. 



Most of the above studies discuss the stationary prop- 
erties of critical systems. On the other hand this knowl- 
edge is often not sufficient for a complete description 
of the dynamics. This need led us in recent years to 
focus on new measures that efficiently probe the com- 
plexity of quantum time evolution. One such measure 
is the so-called Loschmidt Echo (LE) which probes the 
sensitivity of quantum dynamics to external perturba- 
tions (for recent reviews see [l^l)- The recent literature 
on the subject is quite vast and ranges in areas as di- 
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where iJ^ = -\- eW is a one-parameter family of hamil- 
tonians, is the unperturbed hamiltonian, V represents 
a perturbation of strength e and \ipo) is an initial state. 

For a non- critical quantum system with a chaotic clas- 
sical counterpart, the decay of the LE depends on the 
strength of the perturbation parameter e. Three regimes 
have been identified: the standard perturbative, the 
Fermi Golden Rule, and the non-pcrturbative regime. 
The first two can be described by Linear Response The- 
ory leading to a decay which depends on the perturbation 
strength e as F{t) ^ e"''^*-' and F[t) ^ e~'^ *, respec- 
tively [2ll |. In the non-perturbative regime, the LE 
initially follows a Lyapunov decay F{t) ~ e~^*, with a 
rate given by the Lyapunov exponent A of the underlying 



classical system [17|, |18[, whereas for longer times (be- 



yond the so-called Ehrenfest time) the LE decays in ac- 
cordance with the classical autocorrelation function [i^] . 
This behavior matches the decay of the classical echo ac- 
cording to the correspondence principle. 

In this Letter we make the first step in understanding 
the echo decay of dynamical systems with critical chaos. 
The focus of the presentation is on the non-perturbative 
regime where our study revealed a novel result (see Fig. 
1): we have discovered the appearance of an echo gap 
(initial drop of LE) l~ Fg for initial states which are dis- 
tributed in the parts of phase-space where the hamilto- 
nian function exhibits singularities. We have found that 
at the shortest classical time scale (mean free time be- 
tween singular scattering events) Fg scales as 



Fg{a,e,ri) = 



/ci(Xci) where Xci = C • 2— 
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where the sub- indices q/cl indicate the quantum/classical 
scaling fimction for Fg. In Eq. ([2]), a is the order of singu- 
larity of a non-analytical potential, 77 is the characteristic 
spread of the phase space density of the initial classical 
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FIG. 1: Quantum (red line) and classical (black line) LE 
Fit) for the KR with V(q) = log |g| and A'o = 1 in the 
non-perturbative regime: (a) Torus geometry with L = 1, 
A'' — 2^"^ and classical perturbation e = 10~^ (correspond- 
ing to cr ~ 2.0). The width of the initial preparation is ~ 
0.04275. (b) Cylinder geometry, with L = 10^ iV = 2^® and 
classical perturbation e = 0.4 (corresponding to cr ~ 4.17). 
One observes that after the initial Lyapunov decay (inset; 
the dashed line indicates an average Lyapunov decay), F{t) 
follows a power law decay given by the autocorrelation func- 
tion 0, [22]. Here we have rj w 0.15. In both cases we have 
used more than 10^ trajectories for the classical calculation, 
while an averaging over 800 initial Gaussian wavepackets cen- 
tered at different momenta po and go = has been done in 
the quantum calculation. The filled circles indicate the first 
(classical or quantum) nontrivial timestep of Eg = F{n — 2). 
The statistical errors are smaller than the symbol size. 



or quantum state (e.g. its Wigner function) and the con- 
stant C only depends on details of Hq [23|. Moreover, 
we found that the scaling function /(x) behaves asymp- 
totically as /(x ^ 0) ~ X- The above scaling laws were 
derived based on the analysis of classical dynamics and 
confirmed nicely by numerical simulations. 

The apparent deviation of the quantum echo behavior 
from its classical counterpart (see Eq.([21)), can be seen 
as a violation of the quantum-classical correspondence; 
the latter being confirmed in all previous fidelity stud- 
ies. We have found that the origin of this anomalous 
behavior is due to strong diffraction effects which dictate 
the wave dynamics for the class of dynamical systems we 
investigate in this Letter. 

Below we consider a class of parametric Kicked Rotors 
(KR) defined by the time-dependent Hamiltonian 
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where (p, q) is a pair of canonical variables, T and Kq 
are the period and the strength of the kicking potential 
respectively. The class of KRs that we will study below 



FIG. 2: Classical Echo map Eq. (jS]) for the a = 0, /fo = 1 sin- 
gular potential (cylinder geometry with L = 10"^). The initial 
preparation is a box of size 77 = 7r/300 arround (0, 0). The evo- 
lution snapshot is for the shortest nontrivial time scale n — 2. 
The red line is the theoretical prediction of Eq. ((Qj. Upper 
left inset: The Wigner function representation (torus geome- 
try) of the quantum echo map. The light areas correspond to 
negative phase space densities indicating diffraction phenom- 
ena. Lower right inset: A typical phase space for Ko — 1 (20 
trajectories involved up to time 10^ iterations of the map) 



have a potential which is given by 



\q\ for a ^ 0, 
log \q\ for a ~ 0. 
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The a = case corresponds to MIT and will be 
investigated below in detail. The classical dynamics is 
described by the following map: 



Pn+l 



■ p„ - KoV'{qn); qn+l = + Tpn+1, 



(5) 



where all variables are calculated immediately after one 
map iteration and V' = . The domain of q is within 
the interval — tt < q < it. The map ([5]) can be studied 
on a cylinder p d (— c», 00), which can also be closed to 
form a torus of length 27rL, where L is an integer. For 
Ka > the motion is chaotic (see inset of Fig. 2) with 
a (local) Lyapunov exponent given by A(g) = 21og(l -I- 

Ko/i2q') + /fW? + PW(V)P) for q ^ 0. 

The quantum evolution is described by a a one-step 
unitary operator f/g acting on the wavefunction ^/'((j): 

Vq = eyi^{-iTn^ /2) exp(-iA:y(q)) , ?i = 1 (6) 

where n = -id/dq, -N/2 < n < N/2, r = {2tiL/N)T 
and k = [N /2ttL)Kq. Optionally, we define an effective 
Planck constant ?ieff = 2ttL/N. The classical limit cor- 
responds to ^ 00. Without loss of generality we will 
assume below that T = 1. It was shown in that for 
—0.5 < a < 0.5 the eigenfunctions of the above unitary 
operator are multifractal while the levels resemble the 
statistical properties of disordered systems at MIT. 
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FIG. 3: The Fg for various a, rj, e's of the KR defined by 
Eqs. (|3I4|) . In (a) we report the classical Fg against the 
scaling variable Xci (see Eq. ([2}) while in (b) we report the 
corresponding quantum Fg, vs. the scaling variable Xq (see 
Eq. ([2])). At the same sub-figure we report the results for the 
RPKR resulting from the model Eq. (O by a randomization 
of the phases of the kinetic part of Uo- In the insets we report 
a magnification of the main panels in the regime of « 1. 
In all cases an excellent data collapse is observed. 



For the echo calculation, we perturbed our system with 
the following (smooth) potential W{q) = cosg. Cor- 
respondingly the perturbed quantum kicking parame- 
ter is cr = e/hcs- Quantum mechanically, the initial 
preparation is a Gaussian wavepacket centered along 
the line of singularity i.e. {qo,Po) = (0,Po)- The lin- 
ear width 2r] of the packet is taken to be minimal (i.e. 
2r] = Ap = Aq = , where we perform an averag- 

ing over different po's in order to eliminate fluctuations. 
The corresponding classical initial preparation is given 
by a uniform distribution of trajectories located inside a 
box of area A ~ 21] x 2ri ^ hcs. We then define the clas- 
sical LE, Fc\{t ~ n), as the overlap of the initial area Aq 
with the area A f obtained by evolving Aq for n iterations 
of the perturbed map and then reversing the evolution 
for n iterations with the unperturbed one. We have also 
checked that the results remain qualitatively the same 
when we chose an initial classical distribution to be a 
Gaussian density, equivalent to a Wigner function of the 
quantum Gaussian wavepacket. 

We start our analysis with the classical derivation of 
Eq. (U). To this end we consider the classical echo dy- 
namics [26[ . We denote by $o {p, q) the forward sym- 
plectic map defined in Eq. ([5]) while by <i>e = <I>o o 



we denote the corresponding perturbed forward map. In 
this notation P^{p,q) = [p -\- €S:in{q),q) is a symplectic 
map generated by the perturbation (perturbation map). 
In this framework, the n— step echo map is defined as 
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For singular potentials and initial conditions close to 
singularity, we express the phase space shift produced by 
the echo map after the second iteration step: 



Apf = -p^«-esiu(<7i)A'oy"(<7o) 



-<lo ~ -esin(9i) 



(9) 



where qi is given by the map ([5]) with the initial condition 
(Po ' qo) ^ (Poj 9o)- Due to the fact that the initial condi- 
tions are populating a box centered around the singular- 
ity line, we can assume that sin((7i) is a pseudo-random 
variable with density V{x = sin((7)) — (l/7r)(l — x^)^^^^. 
The accuracy of our assumption is tested in Fig. 2, where 
we compare the envelope of the shift in momentum as it 
is given by Eq. ^ with the exact echo dynamics. 

For a given go (and assuming e ^ 77), the probability 
to return to the initial phase- space box is estimated as 



1 



V{qQ) = - 



7rfj + eKQ\V"{qQ)\- 



(10) 



provided that the typical echo shift Ap2 is much larger 
than rj. The echo probability is just the integral of P(go) 
over the initial interval go G [^?7, v]'- 



(11) 



For the specific family of singular potentials discussed 
in this Letter, the above relation gives us: 



^ci - S 1 



2^^010(0-1)1(3-0) 
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for a 7^ 
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Our results in Eq. are nicely confirmed in Fig. 

3, where we are plotting the echo gap for various a, e 
and rj values by making use of the rescaled variable Xci 
given by Eq. ([2]) . Strictly speaking, the above results are 
applicable only for the case where F^y <C 1. Nevertheless, 
our numerics indicates that the scaling behavior Eq. ^ 
continues to apply for values of F^^ ~ 1. 

We have also tested the results of the classical analysis 
against the quantum echo gap. A complete breakdown of 
the quantum-classical correspondence (QCC) is observed 
after the shortest non-trivial time scale (two iteration 
steps) . This is associated with the fact that the Ehrenfest 
time for our system Ie ~ log(/ioff)/A(go) — > when go — > 
[231 . Weak correspondence is restored for longer times 
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when the echo dynamics spreads ergodically resulting in 
a vanishing measure of the critical line at go = 0- 

In Fig. 3 we observe that although the e dependence of 
the quantum Fg is captured by the classical calculations, 
both the ry ~ \/h and the a dependence differ drastically. 
The latter can be explained by a Random Phase Kicked 
Rotor (RPKR) with singular potential, which is simply 
given by Eq. ^ and replacing the eigenvalues of ttt? /2 
by random phases. This indicates that the appearance 
of a gap is insensitive to classical dynamics and thus can 
be captured by a Random Matrix Theory (RMT) model- 
ing which preserves the power-law band structure of the 
evolution operator. Thus the fidelity gap is a universal 
phenomenon of critical systems described by these RMT 
models, and can be used as an alternative criterium to 
level or wavefunction statistics [l|, 0, H, @, 0] ■ 

We argue that the violation of the QCC is due to dom- 
inating diffraction effects that appear as a consequence of 
the singular potential. This is illustrated in the inset of 
Figure 2 where we show the Wigner function (computed 
according to Ref. 28[ ) of the echo map for the first non- 
trivial time step. One observes the appearance of non- 
classical regions in the phase space where the Wigner 
function takes negative values. 

In conclusion, wc find that the LE of dynamical sys- 
tems exhibiting critical chaos decays instantaneously 
with a gap that scales inverse proportionally to the 
strength of the perturbation. The order of the poten- 
tial singularity is encoded in the scaling properties of the 
classical echo gap, while the corresponding quantum gap 
is insensitive to it and its scaling properties are described 
by an RMT modeling of critical systems. This deviation 
is explained in the basis of strong diffraction which is a 
dominating mechanism of short time echo dynamics. 
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